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Abstract: The local polynomial fitting method is applied to estimate the state-price density (SPD) im- 
plicit in financial asset prices. It is shown that the estimator of SPD is consistent under mild 
conditions. Then the VaR (Value at Risk) based on SPD is proposed, which incorporates 
economic valuations, and thus is more reasonable than the traditional VaR. Simulations 
and test of invalidation are done for Black-Scholes models to assess the performance of the 
proposed method. 
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1 Introduction 


Risk management has become a first-order concern of financial managers. The fundamental 
problem in the risk management is how to measure financial risks. For example, the well- 
known Markowitz’s portfolio theory is based on the risk measure characterized by variances of 
risky assets. Many risk measures in finance have been proposed, among which are variances 
or standard deviations, semi-variances, volatility coefficients, market risk @, ô coefficients, VaR 
etc. See Jorion!!! and Penza and Bansal?! for details. It is worth mentioning that VaR now 
becomes very popular in the risk management. Jorion!!! in a milestone book on VaR, defines 
VaR this way: “VaR summaries the expected maximum loss (or worst loss) over a target horizon 
within a given confidence level.” In Ait-Sahalia and Lo’s words!3], “VaR is an estimate, with 
predefined confidence interval, of how much one can lose from holding a position over a set 
horizon.” Indeed, VaR has become one of the most important measure in the financial risk 
metrics. 

Although measures such as variances and VaR do capture the important features of un- 
certainty, they are measures purely in statistical senses and completely ignore the economic 
valuation. To overcome the drawbacks of traditional financial risk measures, Ait-Sahalia and 
Lol?! proposed a new measure, called E-VaR. Based upon the state-price density implicit in 
financial asset prices, the E-VaR inherits many good properties of the traditional statistical 
VaR and at the same time incorporates economic valuation. The reason is that the SPD has 
two characteristics. One is that it can be derived either in preference-based equilibrium models 
or in the arbitrage-based models of Black-Scholes!4) and Merton!®!, The other is that it covers 
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all economically pertinent information, such as investors’ preferences, endowments, asset price 
dynamics, and market clearing etc!!. 

It is crucial to estimate the SPD appropriately by the asset prices in order to calculate the 
E-VaR. There were many discussions about estimation of the SPD. Hutchinson et alll! gave the 
estimator of the SPD by estimating option prices nonparametrically. Rubinstein!” proposed the 
estimator of the SPD by means of implied binomial tree, in which the risk-neutral probabilities 
{rž} are estimated by minimizing the sum of squared deviations between {7} and a set of 
prior risk-neutral probabilities. Ait-Sahalia and Lo!:8] applied the nonparametric and the semi- 
parametric methods to estimate the SPD. They observed that the SPD is the second derivative 
of the price of a call option, and thus differentiated an option pricing formula twice to obtain 
an estimator of the SPD. 

In this paper, we apply the local polynomial method to directly estimate the SPD. As Fan 
and Gijbels!®) pointed out, the local polynomial fitting is superior to the other commonly used 
kernel estimators such as the Nadaraya-Watson (N-W) estimator and Gasser-Miiller (G-M) 
estimator at least in three aspects. Firstly, the N-W estimator leads to undesirable form of the 
bias, while the G-M estimator has to pay a price in variance when dealing with a random design 
model. Secondly, the local polynomial method adapts to various types of designs, i.e., random 
and fixed designs. Finally, the method has no boundary effects: the bias at the boundary stays 
at the same order as that in the interior, without use of specific boundary kernels. Furthermore, 
what makes the local polynomial method appealing is that it applies the least squares principle 
and the regression function, thus the derivatives of the unknown function can be estimated at 
the same time. Here we employ the local polynomials to directly estimate the second derivative 
of the option prices and then give the estimator of the SPD. The performance of the estimator of 
the SPD is also discussed. Ait-Sahalia and Duarte!" also discussed the estimation of the SPD 
by the local polynomial approach, but they did not give proofs of the asymptotic properties of 
the estimators. 

The remainder of this article is arranged as follows. The backgrounds and definition of the 
SPD are introduced in Section 2. In Section 3, we give the consistent estimator of the SPD 
by the local polynomial approach. In Section 4, E-VaR is defined and compared with S-VaR. 
Simulations and test of invalidation are done in Section 5. All proofs are put in Appendix. 


2 The state-price densities 


Implicit in financial asset prices is the state-price density. It gives the price of a security 
that pays one dollar if a given state of nature arises, and zero otherwise. With the suitable 
assumptions for preferences and endowments, under the condition of market completeness which 
ensure that a representative investor with utility function U exists, Ait-Sahalia and Lol! show 
that the date-t price of a security with a single date-T liquidating payoff (wr) is given by 
U'(w T) 


Uw a) fi:(wr)dwr 


P, = E,[Y(wr)Mer] =-[ w Gop 


= 9 [yup fe (ordur È eT- Erpur), (1) 
0 
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where f? (wr) =e"7-) Mt r fe(wr) is called the state-price density (SPD), Mir = Per the 
marginal rate of substitution, and f;(-) the density of the underlying Sr on which the security 
is written, with dS; = p(Sz,t)S¢dt + o (St, t)S:dZz. 

The statistical VaR (S-VaR) is based on the price density f;(-). So it does not reflect 
economic valuations. But the SPD fř(-) covers all economically pertinent information, such 
as investors’ preferences, endowments, asset price dynamics and market clearing. Therefore, 
the VaR based on the SPD (E-VaR) is more appealing than S-VaR. In other words, it is more 
reasonable to use the SPD fř(-) to measure financial risk than the price density f,(-). The 
SPD fř(-) can be applied in the same way as the price density f;(-) in the risk management, 
and it has intimate relations with the price density f;(-). Note that f: is related to the data 
generating process {S+}, whereas f;* aggregates all economically pertinent information. Aït- 
Sahalia and Lol’! summarized their relation as follows: any two of the following imply the 
third: 1) representative agent’s preferences; 2) asset-price dynamics; 3) the SPD. 


3 Local polynomial estimator of the SPD 


It follows from (1) that the date-t price P, of a security with a single liquidating date-T 
payoff (Sr) is given by 


Py = eo"? Bt(Sp) = ee f Y(Sr) f? (Sr)dSr, (2) 


where 7 = T — t. A European call option, written on a stock with date-t price S;, strike price 
X and maturity time T, has payoff Y(Sr) = max(Sr — X,0). Under the assumptions of Black- 
Scholes!“ and Merton’!, according to (2), the date-t price H of the European call option is as 


follows 


E EIE O ont f max(Sr — X, 0)f*(Sr)dSr 
0 


aca ” (Sr — Xfi (Sr)dsr, (3) 


where rt, is the risk-free interests rate between t and T, and ô+, the dividend yield. It is easily 
seen from (3) that oH =e "77 f*(X). Then f(X) = ener FH This indicates that the SPD 
ff (-) is the second partial derivative of the price H of a European call option scaled by the 
compound interests rate. In order to reduce the dimension of the explanatory variables on which 
the option price H depends, Ait-Sahalia and Lol! made the assumption that the SPD depends 
on St, Ttr and 6;,, through the forward price Fyr = S,e("t.7—5.7)7 and the risk-free rate Ter: 
Thus the number of the independent variables of H reduces to four, i.e., H = H (Fir, X, T, Ttr). 

To estimate fř(-), we only need to estimate oH. We follow Fan and Gijbels’ idea about 
local polynomial modeling. Let Z = (Fir, X,T, Tt}, and Zo = (Fioro, X057; Tio,ro) lf H = 
A(Fi,7,X,7,Tt,r7) is smooth and Z is in the neighborhood of Zo, then H can be approximated 
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by 
OH OH 
H( Fir, X, T, Ttr) ~ H(Z — Fi, — m USR, =, 
(Fi, T,Tt,7) (Zo) + OF. pes | tr — Fto,ro) + ax ee Xo) 
OH OH 
Ər led (r= 7) + Ort, gag wt — Tto,r70) 
1 8H 
t3 OF?2, |z=Zo (Fir — Fioro)? 
1 OH 2 1 8? H i 
+t39xX2 Vers E Xo) i 2 Or2 Lose ~ 70) 
1 07H 2 8? H 
20r, zez 067 T Tore) + BR, ax zez Ptr — Ftor) = Xo) 
3H 
ton 70T PA (Fir D Frosro)(T ~ To) 
O° H 
ya OTt r |Z=Zo (Fer — Fto ro) (Ttr — Tto,ro) 
8H 8H 
toxos lon — Xo)(T — To) + Xr gang A Aae = To,r0) 
3’ H 
* Grote, |z=Zo (T = T0) (rt,r = tosto): (4) 


This suggests that a locally weighted polynomial regression can be employed. Let 


Q(B) = X {Pu — Bo — Bi (Fu, — Fer) — Bo(Xi — X) 


i=l 


—B3(ti — T) — Paltri — Tr) — Bs (Furi — Fir)? 

—Be(Xi — X}? — br(Ti — 7)? — BalTt,mi — Ttr)? 

—pbo(Fu ri — Fer) (Xi — X) — bio (Firini — Ft,r)(Ti — T) 

—Bi1 (Furi — Fer )(Ttiri — Ttr) — Pia(Xi — X) (i — 7) 

—PBi3(Xi = X)(Tti ri a Ttr) = Bia(Ti ~ T)(Tti ri = rer) } Kn(Zi = Z), (5) 
where {P;,}? and {Z; = (Fi,,7;, Xi, Ti, Tti) Fi, are realizations of the option price H and its 


explanatory variables Z = (Fyr, X, T, r), B = (bo, A1,--* , G14)’ are parameters, and Ka(-) is a 
four-variable kernel function with bandwidth h. 


Denote by B; (j = 0,1,---,14), the minimizer of (5). Obviously, an estimator of SH is 
25. Therefore an estimator of the SPD fi(-) is given by fè = 2e"t:rT 8g. For convenience of 
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notations, let 


1 Far — Fir San A Timi IET (Fan,r - FR} S (71 — T) (Thr — Tir) 

1 Fer Fir e Tira Tir (Fera — Fer)? © (T2 — T)(fta,r2 — Ttr) 
X = 

L Fina — Far e Tina Ter (Finn — Fer)? e (ta —T) (Team — Ttr) 


W = diag(Kr(Zı — Z), Kn(Z2 — Z), , Kn(Zn ~ Z)), P = (Pa, Post s Pin Y- 
Then (5) can be re-written as 
Q(B) = (P — XBYW(P — XB). (6) 
Minimizing (6) with respect to 3, the estimator Ê of is given by 
Ê = (X'WX)-1X'WP. (7) 


The following theorem gives the consistency of Ê(-) and f¥(-). 
Theorem 3.1 Under conditions (A1)-(A4) in Appendix, it admits that 


AZ) =A +0( (F3) ) +0088), as (8) 
fe(x) = RA +O (EF) ) +O) as (9) 


hold uniformly for Z € D and X € D; where D € Rí and D; € R! are compact. 


4 Economic value at risk 


It follows from Section 3 that when both (S+, t) and o(S:,t) are constants, i.e., u(S;,t) = p, 
a(S, t) =o, the price process 5; of the underlying asset satisfies 


1 2 
St = Soexp { (u- z0 \t+o% }. (10) 
If we assume that {Z:+} is a Weiner process, then S; is a log-normal process, i.e., 
S log(S:/So) — (u — 407)¢]? 
Si ~ ad ex {eee eee i ee cs BL k a1) 
V 2rto S 20°t 


We know from Aït-Sahalia and Loll that under the hypothesis of Black-Scholes"! and 
Merton!®!, the date-t price H of a call option maturing at date T = t+ 7, with strike price X, 
written on a stock with date-t price S; and dividend yield ô+,- is as follows 


162. E E f max(Sr, 0) ff (Sr)dSr 
0 


S:®(d,) oe e "77 X ®(d2), (12) 


i 
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where tai as 
X + ôter +5 
= ce A A A E 2 20 Me d2 = dı — ovr, 
O/T 


with r;,, being a risk-free interest rate between t and T, and ® being a standard normal distri- 


dı 


bution function. It is shown by Ait-Sahalia and Lol?! and Ziegler!!!) that the SPD corresponding 
to the price process {S+} is also a log-normal process, that is 


OH 
* pass Ttr T 
fi (Sr) =e Ox? = 
2 1 log(Sr/St) — (rer — 5,2 — $07)7)? 
~ VS QrraSr exp { = 2027 }, (13) 
while the p.d.f. of Sr conditioned on S; is 
o 1 flog(S7/St) — (u — 307)r]? 
Z7\Zt ~ f(Sr) = Vas. exp { = 5027 \. (14) 


Let U be a utility function of a representative investor and a = (u — r + ô)/o. We define U 
as follows 


a=1: U(Sr) =log (Sp), a#1: U(Sr)=S°%/(i-a). (15) 


It is easily derived that if a representative investor has utility function (15), then his coefficient 
of relative risk aversion is a constant, i.e., p(S7) = a. Especially, for a = 0, he is a risk neutral 
representative investor; for a > 0, he is a risk aversion representative investor. In this case, the 
SPD fř and the corresponding p.d-f. f; are given respectively by (13) and (14). 

With a confidence level a, VaR. based on f; is given by 


P(log(Sr) — log(S;) > VaR | S:) =1—a. (16) 


(16) tells us that the loss of the representative investor is no more than VaR with probability 
1 —a. For example, if a = 5%, it is about 5 times in 100 trials for the representative investor 
to lose more than VaR. Note that 


P(log(Sr) — log(S;) > VaR | S+) = P(Sr > Stexp(VaR) | S+) 


Il 


T fe(Sr)dSr. 


S: exp(VaR) 


Then VaR can be easily determined by (16). Since log(S7/S;) is distributed as N((u — 
40?)r,0?7), VaR is the a-quantile of distribution N((u— }07)r, o?r) 


a 1 exp{ -ETUE Urmia. (17) 


VaR V2TTO 2077 


This VaR only depends on the price distribution f; of the asset and contains no information 
about the representative investor’s risk aversion. In order for VaR to contain both statistical 
and economic valuations, we use the SPD fý instead of the p.d.f. fe in (16). Then we obtain 
the economic VaR-——E-VaR as follows 


P* (log(Sr) — log(S:) > VaR* | S+) =1—a, (18) 
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where the probability is with respect to the SPD f;*. VaR* is called the economic Value at risk. 
It follows from (13) and (18) that VaR* satisfies 


A 1 eeaeee (19) 


VaR* V20TO 2077 


It is easily shown from (13) and (18) that 
VaR* = VaR + (r — ô) — p. (20) 


As we know that S-VaR is only an attempt to provide a single number summarizing the total 
risk in a portfolio of financial assets. It is obvious from (20) that E-VaR has more components or 
information than S-VaR. Moreover, if u = r — 6, which implies a = 0 in (15), then VaR=VaR*. 
Thus the representative investor is risk neutral. If p > r — 6, which means a > 0 in (15) and 
the representative investor is risk aversion. It yields that E-VaR is smaller than S-VaR. Hence, 
what the investor are really thinking about the risk with respect to E-VaR is much “larger” 
than that to S-VaR. The difference between VaR* and VaR is u — (r — ô), which depends on the 
investor’s utility function through a. It indicates from (20) that E-VaR incorporates economic 
informations including the representative investor’s preferences, endowment, market liquidity. 


5 Simulations and conclusion 


To assess the performance of the local polynomial estimator of the SPD, we carried out a 
number of simulations on the Black and Scholes models 


dS; = uStdt + o Std, (21) 


where yz and o are both constants. 

Similarly to that in Ait-Sahalia and Lo!®!, we take u = 7.95%, o = 10.28%, r = 3%, T = 252, 
which means there are 252 trading days in a year, and 7 = 120, 130,--- , 160, respectively. We 
take initial price Sọ = 1 and {Z,} Weiner process. Then from (21), the asset price sequence 
{51, S2,- , Sr} are generated. Based on {s:}7, the option price sequence {P;,}7 are given 
according to formula (12). 

In model (4), when 7 and r are fixed constants, the option price H depends only on the 
future price F;,, and the strike price X. So there are only six parameters 8o, (1, B2, B5, Be, Bo 
in (5), which greatly reduces computations and also alleviates the curse of dimensionality. Here 
the kernel function is taken as K},(r,y) = kn(x)kn(y) where Kn(:) = ¢K(-/A) is one-variable 
Gaussian kernel function 


with h = 0.0198. 

Based on the local polynomial estimator, we computed the E-VaR with confidence level 
a = 0.05. We did 500 simulations in order to make the estimator robust. The result is listed in 
Table 1. It can be seen that the difference between the estimated E-VaR and the true E-VaR 
(Black-Scholes E-VaR) is less than that between E-VaR and S-VaR. Moreover, we calculated 
the invalidation ratio based on (18) from the asset price sequence {.5), S2,--- , Sr}. They are 
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roughly around @ = 0.05 according to different 7. All of these results indicate that the local 
polynomial estimator of SPD behaves quite good. 


Table 1: Output of simulations for Black and Scholes model 


Days to expiration T 


Appendix We make the following assumptions about the model. 
Condition A1 
(a) The probability density function f of Zo is bounded, bounded away from zero and 


Estimated E-VaR 
-0.1342 


Invalidation ratio 


uniformly continuous on Rê. 

(b) The probability density function f,(u,v) of (Zo, Zn)’ is bounded for all n > 1. 

Condition A2 

(a) The kernel K is bounded with compact support. 

(b) |G;(u) — G;(v)| < Cl|u — v|} for all 0 < j < 5 where G;(u) = uK (u) and uJ = 
(ul? u}, ul? ult) with j = ji + j2 + ja + ja, ji 20. 

Condition A3 

The prices H = H(Z) has bounded and uniformly continuous third partial derivatives on 
R4 and these derivatives are Lipschitz continuous. 

Condition A4 

(a) E\P,|” < œ for some v > 2. 

(b) The conditional density fz,)p,(ulv) of Zo given P, exists and is bounded. 

(c) The conditional density f(z,z,)\(P,,Pi4.)(ulv) of (Zo, Zi) given (Pi, Pi+i) exists and is 
bounded for all J > 1. 

(d) The process {P;, Z;} are strongly mixing with the mixing coefficients a(k) satisfying 


So jHal)} T? < 00, 
j=1 
for 2 < y < v and a > 1 — 2/y, while 
a(k) = sup { |P(AB) — P(A)P(B)|: A€ F, Be Perey. 


here F, = F{(P;,Zi),1<i<t}, F* = F{(P;, %),t+k < i}. 


co 


(e) The mixing coefficients a(k) satisfy )>7~, ¢1(m) < oo, here 


Ries n(n) ( 


r1(n) 


nT» 
h4 inn 


4 1/2 
ri(n) = pee Tn = (ninn( In Inn)!*°) 


)a(rs(n)), Li(n) = ( TE y 


1/2 
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Finally, the bandwidth hn satisfies hn — 0 as (n — oo) and 
ntt? pa 
In n[{(Inn)(In Inn)!+4]2/4 


The proof of Theorem 3.1 Let 


— œ, as n—oco forsome 0<6<1. 


k= (k1, k2, ks, ka), k! = kı! x ka! x k3! x kal, 
|k] = kı thotkstka, Zf = 2252239234, 


Íi I 


2 -5S D > yi ky + k2 + k3 + k4 = j, 


alkih(Z) 
O02," AZo" 808Z OZ," 


Then (4) and (5) can be expressed, respectively, as 


(D*h)(Z) = 


H(Z)~% Ð gD*H(Z)|z=29(2 — Zo), (22) 
O<k<2 

Q(8)= 35 (Pa- YD az- EE) (23) 
i=1 O<k<2 


It can be seen from (23) that satisfies the following equation 


tnj(Z)= $ hl*lbx(Z)Snj+e(Z), 0< jl <2, (24) 
os|RI<2 
where 
lc Z-= 253 
tn.j(Z) = a2 Pa -2 Y KAZ - Z), (25) 
1%“ /Zi-Zyi 
Saj = = D> (=) Knl% - 2), (26) 


= 


and K,(u) = #rK(#). Let m = (Tho Tn,1 Tn,2)15x1» Where Tn,o = tno(Z). The first element 
of Tn,ı is the value of tn,;(Z) when j = (0,0,0,1), the second element of Tn,ı is the value of 
tn,j(Z) when j = (0,0,1,0), and so on; the first element of Tn,2 is the value of t,,;(Z) when 
j = (0,0,0,2), the second element of Tn,2 is the value of t,,;(Z) when j = (0,0,1,1), and so 
on. The last element of 7,2 (the tenth element) is the value of t,,;(Z) when j = (2,0,0,0). 
Note that 7,0, 7,1, Tn,2 are counted with respect to |j| = 0,1,2, and tn,; are alphabetically 
arranged according to each j. Similarly, h!*l6,, 0 < k < 2 can be written as an array Gy = 
(Bn,o; Bris Bn,2)! . Note that Bn and B have the same elements except the order. Moreover, all 
possible values of Sn,j+x can be written alphabetically as matrixes S,, |;1 ,;. Let 


Sn,0,0 Sn,0,1 Sn,0,2 
Sn = | Snijo Sni Sn,1,2 


S, S. S, 
n,2,0 n,2,1 n,2,2 15x15 
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Then the equation (24) can be expressed as 
Ta(Z) = Sn(Z)Ba(Z). 
It is easily shown that S» is positive-definite. Thus 


Bn(Z) = Sz (Zyma (Z2). 


Let 
1“ -Zi 
th ;(Z)= sa Pa, — HZ ( f 4 ) Ka(Zi - Z). 
Then 
` Oli% Zi- Z\i ; 
tni(2) = GD = DMZ) AE ) KMA- 2) 
1 
= Š Gr D*H(Z)Snj+r(Z) + eng(Z), 
O<SIk|<2 ~ 
where 


3 n 5 bk ° 
mutt) = 2 EES (BSL) yen 


x A [ (D*H)(Z + w(Z — Z)) — (D*H)(Z) | (0 — w)? dw. 
0 


Note D*H = k!by, then the above results and (24) yield 


th (Z) = JO wR [b.(Z) — bx(Z) ]Snj+e(Z) 


O<|k|<2 


1 
-K X gi P H(Z)Sn,j+4(Z) — en j(Z). 
|kj=3 ` 


367 


(27) 


(28) 


(29) 


Thus, similarly to that for mn, t} ;(Z) can be alphabetically written as following, denoted as 


™(Z) 
TA (Z) = Sn(Z)(B — B)(Z) — h?B,(Z)H3(Z) — en(Z), 


(30) 


where Bn(Z) = (Sn,0,3, Sh,1,3) Sn,2,3), while H3(Z) and en are the matrix and the array, which 
alphabetically written with respect to 3,(D! |H)(Z), |j| = 3 and en, 0 < |j| < 2. Following 


(30), we have 


Bn(Z) — BIZ) = Sr (Zyra (Z) + WPS, *(Z)Bn(Z)H3(Z) + 8, °(Z)en(Z). 


When all above conditions satisfied, it yields from Masry"?! that 


sup | 5ng(Z) -A2 =O (Fa) ) +O) as 


368 CHINESE JOURNAL OF ENGINEERING MATHEMATICS VOL. 27 


where uj = fga vu? K(u)du and 


sup |en ;(Z)| = 0(h3), a.s. 
ZED 


It is shown that (8) holds, and hence (9) holds too. 
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KROSS E EN SBS Wt RHE VaR PHA 
= 7, RARP, PR, fe ee 
(1- NAKA Si BEES. FON 510006; 2- HKZ, MÆ 512005) 


i E: ASC ABS IAA TER AGT oe BE OAR PR RAS OM EE (SPD), ERIA RE 
FHE T SPD ANT AAAI. AEF SPD RIEN SARA OMAN VaR. MiXP EME 
FSH VaR BAH. Sett HM Black-Scholes EHET T BATRA AIRY ETA TPE NES 
KR: RENERE: MSTA: ERME 


